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Abstract 

A  new  type  of  optical  bistability  mediated  by  photon  and  virtual  exciton 
coupling  is  suggested  for  polymers.  Numerical  results  with  parameters 
pertaining  to  polydiacetylene- toluene -sulfonate  are  given  as  an  illustration. 
The  stability  of  the  steady-state  solution  is  examined,  and  the  physical 
mechanism  responsible  for  the  bistability  is  discussed. 
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In  the  transparent  region  well  below  the  absorption  edge  of 

semiconductor  quantum  wells,  a  considerable  amount  of  work  has  been  done  on 

1-4 

nonlinear  coherent  optical  processes  due  to  the  excitation  of  light.  Such 

5-8 

investigations  have  been  extended  to  polymers  in  recent  years.  When  the 

frequency  of  the  incident  laser  beam  is  tuned  below  the  exciton  resonance  in 

the  semiconductor ,  it  is  known  that  virtual  exeifcns  are  generated.  These 

excitons  interact  nonlinearly  with  photons  in  the  same  way  as  real 
2  4-6 

excitons.  ’  This  nonlinear  interaction,  sometimes  known  as  the  phase  space 

filling  effect,  determines  the  spectral  changes  associated  with  the  optical 

excitation  of  two-dimensional  (2D)  systems  of  semiconductor  quantum  wells  or 

of  one -dimensional  (ID)  systems  of  polymers.  It  also  predicts  a  reduction  in 

exciton  oscillator  strength  or  a  bleaching  of  the  exciton  transition. 

6-8 

Several  recent  studies  have  demonstrated  that  it  is  the  optical 

phonon  modes  that  mediate  the  exciton  bleaching  in  nolydiacetylene -  toluene 

sulfonate  (PTS) .  Conjugated  polymers  have  attracted  much  attention  because  of 

(3) 

their  large  nonlinear  optical  susceptibility  \  ■  All  these  point  to  their 

tremendous  potential  of  applications  in  optical  signal  processing  and  high¬ 
speed  computing.  It  is  therefore  of  great  interest  to  investigate  the 
fundamental  me  V-nsra  that  yields  such  a  giant  optical  nonlinearity. 

We  consic.  n  this  Communication  the  special  case  in  which  a  sample  of 
polymer  is  contained  in  an  optical  cavity,  and  study  its  response  to  a  laser 
beam  directed  into  the  cavity.  The  polymer  system  in  our  treatment  is 
described  by  a  Hamiltonian  which  is  obtained  by  adding  the  photon- exciton 
interaction  to  the  model  proposed  in  Ref.  6.  This  simple  model  has  been 
proven  successful  in  the  explanation  of  phonon-mediated  optical  Stark  effect 
observed  experimentally.7  The  purpose  of  this  paper  is  to  investigate  the 
relation  between  the  input  optical  field  intensity  and  the  output  intensity 
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with  special  interest  in  the  possible  existence  of  any  optical  multistability 
in  the  system.  Our  approach  is  semiclassical  in  the  sense  that  Quantum 
fluctuations  are  neglected  as  usual.  The  phonon  and  exciton  variables  are 
then  eliminated  adiabatically  in  our  search  for  the  multistability . 

The  total  Hamiltonian  can  be  written  as 

H  —  fia^a  +  w  b^b  +  co  c^c  +  Ac^c(b  +  b^) 
ph  ex 

+  ig(a^c  -  c^a)  +  i(a^£e  -  aE*e^C)  ,  (1) 

where  we  have  ignored  the  spatial  variation  of  the  cavity  field  which  creates 
the  virtual  exciton.  We  consider  in  (1)  only  one  of  the  many  phonon  modes 
coupling  most  strongly  to  the  exciton  and  neglect  the  momentum  dependence  of 
the  exciton.6  Here,  a^(a),  b^(b)  andc^(c)  are  the  creation  (annihilation) 
operators  for  the  cavity  field,  phonon  and  exciton,  respectively,  with 
corresponding  frequencies  f i,  w  ^  anc^  wex-  ^  is  the  coupling  constant  for  the 
phonon-exciton  interaction,  while  g  stands  for  the  exciton- field  coupling 
constant.  The  coherent  driving  field  is  taken  to  have  the  same  frequency  Q, 
and  its  amplitude  is  E. 

As  mentioned  above,  our  approach  is  semiclassical.  The  procedure  is  as 
follows.  We  first  find  from  the  Hamiltonian  (1)  equations  of  motion  for  the 
operators  a,  b,  c  in  the  rotating  frame.  Then  we  replace  the  operators  by 
their  mean  values  defined  as  classical  variables  a  -  <a>,  -  <b>,  r\  -  <c>  and 
nc  -  <c^c> .  In  other  words,  in  the  limit  of  zero  fluctuations,  all 
correlation  functions  may  factorize  and  we  have  the  optical  Bloch-like 
equations 
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a  **  E  -  ica  +  grj 


(2a) 


*  "  '^ph  '  "  LXnc 


(2b) 


rj  -  -(iA  +  7ex)>7  -  i *1(0  +  0*)  *  go 


(2c) 


nc  -  *g  (<**»?  +  9*a)  -  7exnc 


(2d) 


where  we  have  defined  the  detuning  A  -  w  -  Cl.  The  damning  rates  jc,  ,  ?^d 

ex  ’  pn 

7^^  have  been  introduced  phenomenologically. 

In  the  limit  of  high  Q  cavity,  we  have  k  «  7ex>7ph-  Consequently,  the 

Q 

excitonic  and  phononic  variables  can  be  eliminated  adiabatically .  In  other 
words,  the  time  variation  of  these  variables  are  assumed  to  follow  that  of  the 
fields.  By  setting  the  left-hand  side  equal  to  zero,  it  is  therefore 
straightforward  to  obtain  from  Eqs.  (2b) -(2d) 


0  -  -MWph  +  i7ph)nc/(J;h  +  7ph) 


(3a) 


0  +  0*  =  -2An  /w 

c'  ph 

9  -  -g(l  *  i7)o/(l  +  Y2) 

Apnc  -  2ApAnc  +  (1  +  A2)nc  ’  2§2|q|2  - 

where  we  have  defined 


(3b) 


(3c) 


(3d) 


(4a) 
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A  -  2A2/w  , 
P  Ph 


(4b) 


6,7 


and  have  made  use  of  the  fact  that  u  ^  »  7  ’  All  quantitites  with  the 
dimension  t  ^  are  measured  in  the  unit  7  =  1.  Equation  (3c)  follows 
directly  from  (2c)  and  3(b),  and  (3d)  is  obtained  by  plugging  (3c)  in  (2d). 

For  convenience,  we  now  introduce  the  dimensionless  variables  r  -  ict, 
E  -  E/k  and  a  -  J 2g  a.  The  equation  of  motion  (2a)  for  the  field  mode  then 
becomes 


f7“c  '  y2®2  E  -  ^lOcI2’ 


(5) 


where  the  function  f  is  defined  as 


f(K  I  > 


g  1  -  i[A  -  A  n  (|a  |2)] 

1  +  (— )  - ^ 

1  +  (A  -  Vc(|aJ  )] 


(6) 


That  the  number  of  virtual  excitons  n  is  a  function  of  the  field  intensity  I 

c 

2 

=»  \c^c\  only  is  clearly  seen  from  (3d).  In  the  steady  state,  we  set  ac  “  as  ** 
constant  in  (5)  and  (6).  A  straightforward  calculation  then  leads  to 


1  +  Y  (I .  ) 
_ s  m 


(1  +  g2/«)2  +  Y2(I.n) 


I. 

m 


(7) 


where  we  have  defined  the  steady-state  cavity  field  intensity  ”  I  as  I  •  tbe 

2  2 

driving  field  intensity  I.  -  2g  I E I  and  Y  (I.  )  -  A  -  A  n  (I.  ).  The  fact 
°  in  1  s  m  p  s  in 

that  ng  is  a  function  of  the  driving  field  intensity  I ^  is  seen  from 
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n3  -  2 (A/A  )n2  +  [A2  +  (1  +  g2/*)2]A'2n  -  I.  /A2  -  0  ,  (8) 
s  p  s  6  '  p  s  in'  p  v  1 

which  follows  from  (3d)  and  (7).  Equation  (8)  may  have  one  or  three  real 
roots  depending  on  the  parameters  chosen.  If  there  are  three  possible  values 
of  ns  for  a  given  1^,  Eq.  (7)  implies  that  there  may  exist  three  possible 
intensities  of  the  output  light  since  the  output  intensity  is  proportional  to 
the  cavity  field  intensity. 

It  is  necessary  to  determine  the  stability  of  the  steady-state 
solutions.  This  may  be  done  in  a  similar  fashion  as  discussed  in  Ref.  11. 
Consider  a  general  differential  equation  of  the  type 

®c  "  F(acac)  •  <9> 


which  may  be  linearized  by  a  -standard  procedure  for  small  perturbations  6 a 

c 

and  5a*  about  their  steady-state  solutions  a  and  a*.  The  resulting  linear 
c  s  s  ° 

equation  is 


(6a  I 

'ha  1 

C  -  A 

c 

6a*J 

v  CJ 

5  a* 

L  cJ 

(10) 


where  the  matrix  A  is  given  by 


A 


da 

c 


F  (a 


a*) 

s 


da 

c 


da*  F(as'°*s) 
c 


7^7  F*(a  ,a*) 
5a*  s  s 
c 


(ID 


F*(a  a*) 
s  s 


The  condition  for  stability  is  that  the  eigenvalues  of  A  must  have  negative 
real  parts,  which  requires  that  Tr  j  A |  <  0  and  det|A|  >  0.  More  explicitly,  we 
have 


k (ReF>  -  k  (ImF)  < 0 


(12a) 


k  <ReP)  k  (I"F)  ■  k  (EeF)  k  (I"F)  >  0 


(12b) 


where  u  -  Re(a  )  and  v  -  Im(a  )  .  Thus  the  stability  conditions  for  our  Eq 
(5)  are 


f  1  (I  )  +■  I  77  f,  (I  )  >  0 
Is  soils 


I 


1-J. 


f  (I  )[f  (l_)  +  21  rr  f  (I  )]  >  0 

IS  1  5  Solis 


(13a) 


(13b) 


where  £^(£^)  is  the  real  (imaginary)  part  of  f. 

On  the  other  hand,  since  Eq .  (8)  has  oniy  one  real  root  for  A  <  0,  it  is 

quite  clear  that  there  is  no  multiple  output  intensity  whenever  the  cavity 

field  is  tuned  on  or  above  the  resonance  with  the  exciton,  no  matter  how 

strong  the  input  intensity  is.  The  situation  is  tne  same  wnen  no  phonons 

are  involved  or  when  A  =-0.  Therefore,  the  nresence  of  the  optical 

P 

multi  stability  is  really  mediated  by  the  phonons  via  their  nonlinear 
interaction  with  virtual  excitons . 

It  is  known  that  polydiacetylenes  possess  not  only  large  but  also 
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small  transmission  loss  a. 


An  added  attraction  is  their  flexibility  in  the 


We  consider  in  what  follows  polydi acetylene  - 


construction  of  waveguides. 

to  luene  -  sul  f  onate  (PTS)  as  a  specific  example  in  our  numerical  calculation. 

Following  Ref.  6,  we  take  only  one  of  the  phonon  modes  which  couple  most 

strongly  to  the  excitons.  The  parameters  for  PTS  are^’^  7^  -  0.05  eV,  w  ^  - 

0  2^8  eV,  A  -  0.1  eV  and  7  ,  -  0.002  eV.  The  steady-state  intensity  I  has 

ph  s 

been  calculated  numerically  as  a  function  of  I.  with  the  energy  unit  7 

m  ex 

2 

Thus  «  .  -  5.]  6  and  A  -  2.0.  The  results  fcr  different  values  of  A  a.nA  g  /x 
ph 

are  plotted  in  Fig.  1.  It  is  found  that  only  when  the  cavity  field  is  tuned 

sufficiently  below  the  exciton  resonance,  say  bv  an  amount  >  ~  w  .  ,  does  there 

'  '  -  l  ph 

exist  optical  tristability .  'When  the  cavity  damping  is  large  or  when  g“/x  is 

small,  there  is  no  such  nonlinear  phenomenon  either.  This  means  that  the 

cavity  feedback  is  one  of  the  vital  factors  to  creat  optical  tristability  in 

the  present  case.  On  the  other  hand,  such  nonlinear  phenomenon  does  net 

appear  in  cavities  of  very  high  quality,  namely,  cavities  of  exceedingly  high 

2 

figure  of  merit  Q,  or  g  /x  is  large.  This  implies  that  enough  line  width  of 
the  cavity  is  required  to  create  such  nonlinear  interactions,  and  hence  that 
this  tristability  is  at  least  a  dispersive  one. 

Through  the  linearized  stability  analysis  of  the  steady-state  solutions 
discussed  above,  we  have  found  that  all  the  multiple  solutions  are  stable 
except  for  those  in  the  vicinity  of  turning  points  such  as  A  ar.d  B  in  Fig.  1 


(b)  and  (c).  Hence,  when  the  input  field  intensity  increases  continuously, 


passing  the  turning  point,  and  then  decreases,  the  output  field  intensity  will 

follow  a  reversed  hysteresis  loop.  In  contrast  to  the  ordinary  optical 

bistability  behavior,  this  novel  multistability  behaves  more  like  the  special 

hysteresis  loop  in  systems  that  exhibit  increasing  absorption  optical 

13 

bistability.  It  is  also  of  interest  to  note  that  tne  middle  branch  as  well 

as  the  other  two  as  shown  in  the  figure  is  stable.  This  is  of  course  very 
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different  from  what  happens  in  ordinary  cases.  However,  as  long  as  the  input 

field  intensitv  .ranges  monotoni«.aily ,  the  cavity  field  intensity  is  not 

expected  t-  „ove  into  the  middle  branch.  For  instance,  the  system  follows  the 

upper  branch  as  I  increases  until  it  reaches  the  turning  point  A  where  it  is 

unstable  and  jumps  to  the  lower  branch  as  I.  increases  further  Similarly 

m 

the  svstem  follows  the  lower  branch  and  jumps  to  the  upper  branch  at  the 
unstable  point  B  when  I  decreases.  It  is  not  impossible  to  reach  the  middle 
branch  if  the  input  intensity  reverses  its  course  of  change  at  the  turning 
points.  In  addition,  we  learn  from  the  numerical  studv  that  the  threshold  for 
" he  bistability  is  raised  by  increasing  the  detuning  A.  Thus,  a  proper  choice 
of  A  may  be  preferable  for  the  realization  of  optical  bistability. 


r«e  now  attempt  to  understand  physically  the  mechanism  responsible  for 
this  novel  phenomenon  of  reversed  optical  bistability.  It  is  seen  from  Eq . 
(3c)  that  the  resonance  enhancement  of  the  dipole  moment  in  steady  state 


occurs  at  u  -  fi  - 
x 

relaxation  energy. 


An  =0.  The  term  A  n  mav  be  regarded  as  the 
p  s  p  s  ° 

Therefore  both  the  real  and  imaginary  parts 


lattice 
of  the 


optical  vesponse  of  the  PTS  system  depend  upon  this  relaxation  energy,  which 


is  directly  proportional  to  the  number  of  virtual  excitons  induced  bv  the 


cavity  field.  Since  both  the  absorption  and  dispersion  of  PTS  are  known  to  be 

intensity  dependent,  it  appears  that  the  bistability  discovered  in  the  present 

work  belongs  to  a  mixed  type  of  dispersion  and  absorption.  As  we  have  found 

in  our  numerical  analysis,  when  I.  increases  from  zero  and  reaches  a  point 

in  r 

where  n  can  have  three  different  values,  the  largest  n  always  satisfies  the 
s  6  s 

relation  A  -  A  n  -  0.  Further  increase  of  I.  beyond  that  point  leads  to  an 
p  s  m  '  r 

increase  of  the  refractive  index  of  ^TS  .  This  may  be  understood  in  the 
following  way.  As  the  input  intensity  increases  beyond  a  certain  level,  the 
number  of  virtual  excitons  becomes  large  enough  that  the  lattice  relaxation  is 
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sufficiently  high  for  the  refractive  index  to  deviate  from  the  cavity 
resonance.  Consequently ,  the  cavity  field  switches  to  the  lower  branch.  The 
dispersion  effect  plays  a  crucial  role  in  this  type  of  optical  bistability. 

As  concluding  remarks,  we  note  that  the  virtual  excitons  and  lattice 
vibrations  in  PTS  indeed  damp  much  faster  than  the  damping  in  an  ordinary 
cavity.  This  enables  us  to  eliminate  the  excitonic  and  phononic  variables 
adiabatical ly  for  the  special  kind  of  bistability  discussed  above. 
Furthermore,  the  novel  phenomenonon  of  optical  tristabilitv  is  found  by 
considering  only  a  single  phonon  mode  in  PTS.  It  is  therefore  expected  to 
find  optical  multistability  if  more  phonon  modes  are  included. 
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Figure  Caption 

1.  Intensity  of  the  cavitv  field  I  versus  the  input  field  intensity  I.  . 

s  r  7  in 

(a)  A  -  -  5.16,  g2/*  -  0.1;  (b)  A  -  ,  g2/ic  -  0.5;  (c)  A  -  w 

g2/*  "  1-0;  (d)  A  -  g2/*  -  1-5;  (e)  A  -  wph>  g2/*  -  2.0;  (f)  A  -  3.0, 

2 

g  / k  -  1.0.  The  dashed  lines  in  (b)  and  (c)  are  hypothetical  to  complete 
the  reversed  hysteresis  loop. 
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